Arbitrary accuracy iterative phase estimation algorithm as a two qubit benchmark 
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We discuss the implementation of an iterative quantum phase estimation algorithm, with a single 
ancillary qubit. We suggest using this algorithm as a benchmark for multi-qubit implementations. 
Furthermore we describe in detail the smallest possible realization, using only two qubits, and 
exemplify with a superconducting circuit. We discuss the robustness of the algorithm in the presence 
of gate errors, and show that 7 bits of precision is obtainable, even with very limited gate accuracies. 



Solid-state quantum computing is now entering the 
stage of exploration of multi-qubit circuits. Co- 
herent two-qubit coupling has been experimentally 
realized for all major types of superconducting 
qubits^'&^'&^ i 8 ' 9 ' 1 ' 1 1 , and two-qubit gates have been 
demonstrated for charge 8 -, phas o 9 ' 10 and flux qubitsii. 
The question then arises, what kind of testbed appli- 
cation can be performed having at hand a very limited 
amount of qubits? 

Here we propose to employ the Phase Estimation Algo- 
rithm (PEA), which can be implemented with just two 
qubits. Furthermore, we suggest how to use this algo- 
rithm to characterize (benchmark) qubit circuits. The 
PEA is an algorithm to determine the eigenvalue of a 
unitary operator U; it is closely related to the Quantum 
Fourier Transform (QFT), which is a key element of many 
quantum algorithms, e.g., Shor's factoring algorithm 12 
and in general Abelian Stabilizer type of problems^. The 
algorithm's relevance for quantum simulations was no- 
ticed by Abrams and LloycU^, and recently emphasized 
by Aspuru-Guzik et ai^ simulating quantum compu- 
tation of the lowest energy eigenvalue of several small 
molecules. It is clear that the PEA will be one of the im- 
portant algorithms in future quantum information pro- 
cessing applications, and how accurately a phase can be 
determined will be an important figure of merit for any 
implementation. 

The textbook 1 - implementation of this algorithm re- 
quires n qubits representing the physical system in which 
U operates, and m ancillary qubits for the work regis- 
ter. The number m determines the algorithm's precision 
1/2™, i.e the number of accurate binary digits extracted. 

There is also an alternative algorithm proposed by 
Kitaev-i2,, where the Fourier transform is replaced with 
a Hadamard transform. In implementing this algorithm 
to obtain a precision of order 1 /2 m , it is possible to run ei- 
ther log m rounds (iterations) with m ancillary qubits or 
m log(m) rounds with only a single ancilla. The precision 
increases exponentially with the number of rounds, but 
each round requires exponentially many applications of 
U, unless powers U 2 are available by different meansil. 

Also the QFT-based PEA can be implemented in a 
multiround fashion, using a single ancillary qubit, based 
on the semiclassical QFT 18 . In this paper, we refer to this 
single ancilla QFT based PEA as iterative PEA (IPEA). 
The iterative version of Kitaev's algorithm is referred to 



as Kitaev's PEA. The relevance of the iterative PEA as a 
viable alternative to the textbook version was noticed by 
Mosca & Ekerfei 9 - in the context of the hidden subgroup 
problem, by Zalka32, for factoring, and by Childs et. al^- 
and Knill et. al. 22 in more physical contexts. 

As long as the number of qubits is a limiting factor, 
implementations of phase estimation with only a single 
ancillary qubit will be of foremost importance. Thus, it 
is instructive to compare the iterative PEA with Kitaev's 
PEA. In the IPEA scheme, the bits of the phase are 
measured directly, without any need for classical post- 
processing. Moreover, each bit has to be measured only 
once, compared to log(m) times. When the phase has 
a binary expansion with no more than m bits, the IPEA 
detcrministically extracts all bits, in contrast to Kitaev's 
PEA which is always probabilistic. The IPEA is also op- 
timal in the sense that a full bit of information is gained 
in each measurement^. 

Theoretically the accuracy of the algorithm is limited 
only by the number of rounds, but in practice it will be 
limited by experimental imperfections. Thus, the exper- 
imentally maximally obtainable accuracy can serve as a 
benchmark for any multi-qubit implementation. 

For benchmarking purposes a setup is needed where 
the phase to be measured can be set to an arbitrary value. 
We describe in detail such an implementation in a system 
of two superconducting qubits. Introducing gate noise, 
we also perform a robustness analysis, indicating which 
gates are most critical, and we calculate the number of 
repetitions needed as a function of noise levels. 

Iterative PEA. We now describe the IPEA briefly, but 
still in some detail, in order to make the robustness anal- 
ysis clear. The most straightforward approach for phase 
estimation is shown in Fig. [TJ The upper line is the ancil- 
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FIG. 1: Naive implementation of the phase estimation algo- 
rithm. 

lary qubit which is measured, and the lower line describes 
the qubits representing the physical system in which U 
operates. Initially the ancillary qubit is set to |0) and the 
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lower line register to an eigenstate |\&) of the operator U 
with eigenvalue e t27 "^ . Right before the measurement the 



system state is \ [(l 
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giving the probability Pq — cos 2 (7T0) to measure " 0" . 
By repeating this procedure N times, Pq can be deter- 
mined to an accuracy of 1/y/N. Thus, one needs to go 
through at least N ~ 2 2m independent rounds to obtain 
to accurate binary digits of 0. The number of rounds 
corresponds to the number of measurements since each 
round is terminated with a measurement. 

Kitaev's PEA allows the number of rounds and conse- 
quently the number of measurements to be drastically re- 
duced, with the assumption that the controlled-C/ 2 gates 
are available^. For each k, 1 < k < m, the controlled- 
U 2 gate is used to prepare an ancillary qubit in the 
state ^= (JO) + e*^" 1 *) |1)) . After a number of repe- 
titions, the ratio of resulting zeros and ones is used as an 
estimate for the fractional part of 2 fc_1 0. A classical algo- 
rithm with polynomial runtime is then used to assemble 
from the fractional parts. The whole algorithm performs 
estimation of with precision l/2 m+2 and error proba- 
bility < e after 0(m\og(m/e)) measurements. The gate 
U is applied 0{2 m log(m/e)) times to create the powers 
U 2 , which is nearly a quadratic improvement compared 
to the naive version of the PEA. 

The iterative PEA differs by the following modification 
of the above described procedure: first less significant 
digits are evaluated and then the obtained information 
improves the quantum part of the search for more sig- 
nificant digits. The information transfer is done with an 
extra single qubit Z-rotation that is inserted into the cir- 
cuit, as shown in Fig. [2] Note that k is iterated backwards 
from m to 1. 
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FIG. 2: The fcth iteration of the iterative phase estimation 
algorithm. The feedback angle depends on the previously 
measured bits through LUk = —2Tr(0.0xk+iXk+2---Xm), and 

= 0. 

In order to derive the success probability for each bit 
being determined correctly, we first assume the phase 
to have a binary expansion with no more than m 
bits, — (0.0102 ■■■ 0mOOO .. .). In the first iteration 
(k = m) a controlled-C/ 2 ™ 1 gate is applied, and the mth 
bit of the expansion is measured. The probability to 
measure "0" is Pq = cos 2 [7r(O.0 m OO ...)], which is unity 
for m = and zero for m = 1. Thus, the first bit 
TO is extracted deterministically. In the second itera- 
tion (k = ni — 1) the measurement is performed on the 
(to — l)th bit. The phase of the first qubit before the 
Z-rotation is 27r(O.0 m _i0 m OO . . .), and performing a Z- 
rotation with angle Lo m -\ = — 27r(O.O0 m ), the measure- 



ment probability becomes Pq = cos 2 [7r(O.0 m _iOO . . .)]. 
Thus, using feedback the second bit is also measured de- 
terministically, and generally using the feedback angle 
u>h = —27r(0.0(j)k+i4>k+2 ■ ■ ■ 0m) all to bits of are ex- 
tracted deterministically. 

Denoting the first to bits of the binary expansion of 
the phase as — 0.0102 ■ ■ • 0m, there is in general a 
remainder < 6 < 1, defined by = + <52~"\ In this 
case, the probability to measure m is cos 2 (nS/2). If 
m was measured correctly, the probability to measure 
TO _i in the second iteration is cos 2 (7r£/4) , and so on. 
Thus, the conditional probability Pj~ for each bit to be 
measured correctly is Pk = cos 2 (-7r 2 k ~ m ~ 1 5), and the 
overall probability for the algorithm to extract is 



m = n Pk 



sin 2 (7T S) 



fc=i 



2 2m sin 2 (Tr2- m S) 
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which is the same outcome probability as the textbook 
phase estimation, based on the QFT— . For S < 1/2 
the best m-bit approximation to is indeed 0, while 
for S > 1/2 rounding up to + 2~ m is better. The 
probability to extract + 2~ m is P(l — 6). The success 
probability P(5) decreases monotonically for increasing 
to. In the limit to — > oo, we find the lower bound for the 
probability to extract the best rounded approximation 
to as P(l/2) = 4/-7T 2 . Best rounding implies an error 
smaller than 2~( m+1 ), while an accuracy of 2~ m implies 
that we accept both answers + 2~ m and 0. The success 
probability is then P(<5) + P(l — S), with a lower bound 
of 8/ ir 2 . In conclusion, iterative PEA determines the 
phase with accuracy l/2 m and with an error probability 
e < 1 — 8/tt 2 , which is independent of to. 

Success probability amplification for the textbook 
PEA was discussed in Ref. 13. When the algorithm is 
executed with to' = to + log (2 + l/2e) ancillas, the esti- 
mate is accurate to to bits, with probability at least 1 — e. 
A similar approach can also be used in iterative PEA, by 
determining to' bits and keeping only the to most signif- 
icant bits. However, implementing the U 2 gate for large 
k is the algorithm's bottleneck in a realistically noisy en- 
vironment, as discussed below. 

A different approach, avoiding this problem, is to re- 
peat the algorithm a number of times, choosing the most 
frequent result. In natural ensemble systems, such as 
NMR, this is naturally exploited with advantage^- 4 . How- 
ever, for single systems such as superconducting qubits, 
repetitions on all to bits are unnecessarily expensive. 
From Eq. |T]) it is clear that the main contribution to the 
error probability comes from the least significant bits. 

Thus, we may lower the error probability significantly 
by repeating the measurement of only the first few bits a 
limited number of times, and using simple majority vot- 
ing. It is clear from the binomial distribution that the 
bitwise error probabilities decrease exponentially with 
the number of repetitions. Also, because of the feedback 
procedure the bare error probability sin 2 (7r 2 k ~ m ~ 1 5) al- 
ready decreases exponentially with decreasing k. Thus, 
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one needs only O [log 2 (1/e)] extra measurements to ob- 
tain an error probability smaller than e, independently 
of m. 

Benchmark circuit. The minimal system for imple- 
menting the iterative PEA is a two qubit system, where 
one qubit is a read-out ancilla, and the second qubit rep- 
resents a physical system. From the work of Barenco et 
al^ we know an explicit construction of any controlled- 
U gate, where U is an arbitrary single qubit gate. 
This construct involves three single qubit gates and two 
controlled-NOT (CNOT) gates. 

For benchmarking purposes we propose to use the very 
simple Z-rotation operator 



U 











(2) 



where a is an arbitrary rotation angle. 

The advantages of this operator are 1) it is diago- 
nal in the qubit eigenbasis, thus the initial preparation 
of its eigenstate is straightforward, 2) the phase to be 
measured is controlled directly, and 3) controlled pow- 
ers of this gate are generated by a single entangling 
gate, ZZ(a) = diag (e" iQ , e la , e la , e" M ); this gate can 
be straightforwardly implemented by using most com- 
mon superconducting qubit coupling schemes. As shown 
in Fig. [31 a step of the iterative PEA is implemented us- 
ing one ZZ-gate, and in addition only three single qubit 
gates. The phase we are measuring in this case is set by 
the coupling strength A, rather than by the free qubit en- 
ergy that is the case using the general construction with 
two CNOT gates, a2 k - 1 = AT, T is the pulse duration. 
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FIG. 3: A gate sequence implementing the fc-th step of the 
iterative phase estimation algorithm, on a two qubit system 
using the entangling gate ZZ(a). 

Let us consider implementation of the ZZ-gate with 
superconducting qubits in more detail. For supercon- 
ducting charge and charge-phase qubits operated at the 
charge degeneracy point, and physically connected via 
a Josephson junction placed at the intersection of the 
qubit loop-shaped electrodes, inductive interaction of 
persistent currents circulating in the loops creates direct 
switchable zz-couplin g 26 i 27 i 28 . Thus the implementation 
of the ZZ-gate is straightforward. 

Furthermore, the ZZ-gate is a generic gate for the 
qubits coupled via tunable linear oscillator: this gate 
is generated by applying a composite dc-pulse sweeping 
through the qubit-cavity resonances as shown in 2 -. 

For the permanent transverse coupling (xx-coupling in 
the qubit eigenbasis) frequently discussed in the context 
of charge^, phased, and flux 30 qubits, the ZZ-gate can be 
realized with dynamic control scheme o 30 ' 31 . The para- 
metric coupling method^, suggests harmonic modulation 



of the coupling strength X(t) with the two resonant fre- 
quencies corresponding to the sum and the difference of 
the qubit energies. This induces the Rabi rotations, U^, 
and Uft, in the parallel spin (|00), |11)), and antiparallcl 
spin ( 1 01 } , 1 10)) subspaces, respectively. The ZZ-gate is 
then obtained, up to a single qubit Z-rotation, by apply- 
ing the Hadamard gates to both the qubits, H = Hi C3H2, 
according to the scheme, ZZ(a) = H [U^(a)®U^(a)] H. 

The method 31 is similar although more time consum- 
ing. In this case, the resonant rf-pulses are simultane- 
ously applied to both the qubits, inducing Rabi rotations 
U R . The pulse amplitudes are set equal to the half dif- 
ference between the qubit energies. Such an operation 
produces the gate which is equivalent to the rotation in 
the 1 00), 1 11) subspace, U R (a) = HC/£(a/4)H. Rotation 
in the |01), 1 10) subspace can be performed by first swap- 
ping the states of one of the qubits, and then apply- 
ing the same pulse, E/|(a/4) = X{SU R {a)EX 1 . Thus 
the ZZ-gate is achieved applying the Rabi pulses twice, 
the full operation sequence taking the form, ZZ(a) — 
U R (ia)Z 1 U R (ia)Z 1 . 
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FIG. 4: The success probability of the IPEA to correctly de- 
termine the phase a, with precision better than 2 -5 (green) 
and 2~ 7 (blue) as a function of the noise level. The three cases 
of pure A-gate errors (dotted), pure dephasing (dashed), and 
both types of noise acting simultaneously (full) are consid- 
ered. The simulation was averaged over a evenly distributed 
in the range — n < a < n. 

Robustness analysis. There are numerous imaginable 
sources of error, in all parts of the algorithm from ini- 
tialization via gate manipulation to readout. With our 
setup initialization will probably be accurate, but the 
gates will certainly suffer from imperfections due to en- 
vironmental noise. First we consider the effect of pure 
dephasing with rate T v , which eventually will limit the 
accuracy of the ZZ(a2 k ) gate. In addition, we consider 
imperfect ^-rotations of the form R x (±ir/2 + S x ), where 
5 X is a normally distributed random angle with variance 
A x . These errors modify the probability of measuring 
the correct value of the feth bit into 
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FIG. 5: The total number of measurements needed to obtain 
the phase a with precision better than 2 _m (2 < m < 11), 
with an error probability e < 0.05. 



In Fig.[3]the algorithm's success probability, as a func- 
tion of the dephasing rate and variance is shown for 
m = 5 and m — 7. The algorithm is rather robust against 
x-rotation errors, while being much more sensitive to de- 
phasing, which is evident from the exponentially growing 
factor 2 k a in the exponent of Eq. ([3]). 

As discussed below Eq. (p}, the success probability can 
be improved by repeated measurements of each bit. To 
achieve an overall success probability of 1 — s, we need 
to increase the per bit success probability to (1 — e/m), 
using Nk repeated measurements. For P' k close to 1/2, 



many repetitions are needed, and the binomial distribu- 
tion approaches the normal distribution giving, 



N k = 



1 (er]- x [\ 
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where er/ _1 is the inverse of the error function. Consid- 
ering the dominating effect of dephasing, we find that 
the number of repetitions grows quickly with the de- 
sired number of bits m, Nk oc e 2 ' Q l 2 r v/ A . In Fig. Owe 
plot the total number of measurements N to t — J2k^k, 
needed to obtain 2 < m < 11 bits of the phase a, with 
an error probability e < 0.05. For a realistic dephas- 
ing rate of 1 to 10 percent of the qubit-qubit coupling 
(0.01 < r^/A < 0.1), between 5-8 binary digits of a can 
be extracted with less than 10 4 measurements. 

In conclusion, we have described a benchmark imple- 
mentation of the iterative PEA on two superconducting 
qubits, and analyzed its robustness towards dephasing 
and gate errors. The number of extractable binary dig- 
its is mainly limited by the dephasing, and for realistic 
parameters amounts to 5-8. We believe phase estimation 
will be an essential part of future applications of quan- 
tum computing, and propose that the number of accurate 
binary digits can be used as a benchmark for multi-qubit 
implementations . 

This work is supported by the European Commission 
through the IST-015708 EuroSQIP integrated project 
and by the Swedish Research Council. 



1 J. B. Majer et al, Phys. Rev. Lett. 94, 090501 (2005). 

2 H. Xu et al, Phys. Rev. Lett. 94, 027003 (2005). 

3 A. C. J. ter Haar, PhD thesis, Delft University (2005). 

4 M. Grajcar et al, Phys. Rev. Lett. 96, 047006 (20 06). 

5 S. H. W. van der Ploeg et al, (2006); |cond-mat/0605588"l 

6 A. O. Niskanen, K. Harrabi, F. Yoshihara, Y. Nakamura, 
and J. S. Tsai, Phys. Rev. B 74, 220503(R) (2006). 

7 Yu. A. Pashkin, T. Yamamoto, O. Astafiev, Y. Nakamura, 
D. V. Averin, and J. S. Tsai, Nature 421, 823 (2003). 

8 T. Yamamoto, Yu. Pashkin, O. Astafiev, Y. Nakamura, 
and J. S. Tsai, Nature 425, 941 (2003). 

9 R. McDermott et al, Science 307, 1299 (2005). 

10 M. Steffen et al, Science 313 1423 (206). 

11 J. Plantenberg, P. C. de Groot, C. J. P. M. Harmans, and 
J. E. Mooij, Nature 447, 836 (2007). 

12 P. W. Shor, SIAM J. Comput. 26 1484 (1997). 

13 A. Yu. Kitaev, Electronic Colloquium on Computational 
Complexity 3(3) (1996). 

14 D. S. Abrams and S. Lloyd, Phys. Rev. Lett. 83, 5162 
(1999). 

15 Alan Aspuru-Guzik et al, Science 309, 1704 (2005). 

16 R. Cleve, A. Ekert, C. Macchiavello and M. Mosca, Proc. 
R. Soc. London A, 454, 339-354 (1998); Michael A. Nielsen 
and Isaac L. Chuang, Quantum Computation and Quan- 
tum Information, Cambridge University Press (2001). 

17 V. Giovannetti, S. Lloyd, and L. Maccone, Phys. Rev. Lett. 
96, 010401 (2006). 



18 R. B. Griffiths and C.-S. Niu, Phys. Rev. Lett. 76, 3228 
(1996). 

19 M. Mosca and A. Ekert, Lecture Notes in Computer Sci- 
ence 150 9, 174 (1999). 

20 C. Zalka, |arXiv:quant-ph/9806084"l 

21 A. M. Childs, J. Preskill, J. Renes, J. Mod. Opt. 47, 155- 
176 (2000); 

22 E. Knill, G. Ortiz, R. D. Somma, Phys. Rev. A 75, 012328 
(2007). 

23 W. van Dam, G. M. D'Ariano, A. Ekert, C. Macchiavello, 
and M. Mosca, Phys. Rev. Lett. 98, 090501 (2007). 

24 J. A. Jones and M. Mosca, Phys. Rev. Lett. 83, 1050 
(1999). 

25 A. Barenco et al, Phys. Rev. A, 52, 3457 (1995). 

26 J. Q. You and F. Nori, Phys. Rev. B 68, 064509 (2003). 

27 G. Wendin and V. S. Shumeiko, in Handbook of Theor. and 
Comp. Nanotech., eds. M. Rieth and W. Schommers, ASP, 
Los Angeles, (2006), Vol. 3, pp. 223-309. 

28 J. Lantz, M. Wallquist, V. S. Shumeiko, and G. Wendin, 
Phys. Rev. B 70, 140507(R) (2004). 

29 M. Wallquist, V.S. Shumeiko, and G. Wendin, Phys. Rev. 
B 74, 224506 (2006). 

30 P. Bertet, C. J. P. M. Harmans, and J. E. Mooij, Phys. 
Rev. B 73, 064512 (2006). 

31 C. Rigetti, A. Blais, and M. Devoret, Phys. Rev. Lett. 94, 
240502 (2005). 



